Nonequilibrium phonon backaction on the current noise in atomic-sized junctions 
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We study backaction effects of phonon heating due to tunneling electrons on the current noise in 
atomic-sized junctions. Deriving a generalized kinetic approximation within the extended Keldysh 
Green's functions technique, we demonstrate the existence of a characteristic backaction contribution 
to the noise in case of low external phonon damping. We provide a physically intuitive interpretation 
of this contribution at large voltage in terms of slow fluctuations of the phonon occupation, and show 
that it generally gives a significant correction to the noise above the phonon emission threshold. 
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Introduction. Inelastic transport spectroscopy is an 
important tool of investigation for a wide range of phe- 
nomena in nanojunctions, ranging from vibrations in 
atomic wiresi to magnetic excitations of aggregates on 
surfacesi^ Inelastic studies provide both direct informa- 
tion on the non-electronic degrees of freedom, such as 
vibrations and/or spin excitations and, indirectly, addi- 
tional knowledge of the electronic subsystem. Nowadays, 
experimental interest is not only limited to the mean 
current but it also extends to the current noise which 
provides further information on the transport properties 
of the junction, such as the number and transparency 
of the transmission eigenchannelsi^ Inelastic effects on 
the current noise are also presently under experimen- 
tal investigation;^ accompanied by an intense theoretical 
activityj^"— 

Since the electron-phonon (e-ph) interaction is often 
very weak, a perturbative approach to the lowest-order 
in e-ph coupling is usually appropriate for evaluating 
both the curren1>i2rii and noise ^"—lii enabling semi- 
analytical treatments that can be combined with ab ini- 
tio calculations j^^^ There is, however, a tricky aspect 
of the problem, related to the heating of the vibrational 
mode(s) due to tunneling electrons, whose experimental 
signature is the voltage dependence of the conductance 
above the phonon emission thresholdJ^ Accounting for 
this effect pushes the theory well beyond the plain lowest 
order perturbation expansion. 

Nonequilibrium phonon heating has been taken into 
account either phenomenologically in terms of a rate 
equation for the phonon occupation^iii or microscop- 
ically by infinite resummation of the electron-hole po- 
larization bubbleJ^ii^ As far as only the current is con- 
cerned, the two approaches yield essentially identical re- 
sults in the limit of weak e-ph coupling. However, they 
differ significantly in the case of the current noise due 
to the feedback of the phonon dynamics on the statistics 
of the transmitted electrons, which is not captured by 
the phenomenological rate equation approach. Accord- 
ing to Ref. [13, this feedback results in a cubic increase of 



the noise as a function of voltage at large bias, S ^ . 
However, this prediction contradicts a preliminary un- 
published work by Jouravlev supervised by Nazaro\J^ 
suggesting S ^ V^. A clear physical interpretation of 
the results was not offered in either case. 

Thus, despite its relevance for several systems of ex- 
perimental interest ^li^i^S the effect of phonon heating 
on the current noise still remains an open problem. In 
this work, we develop a consistent microscopic theory to 
address this issue. We derive an analytic kinetic-like ex- 
pression for the phonon Green's function (GF), which 
encodes the effects of electronic fluctuations and allows 
one to take into account the backaction of nonequilib- 
rium phonons on all current cumulants. Focusing on the 
current noise, we demonstrate that current-driven fluc- 
tuations of the phonon occupation lead to a distinct cor- 
rection to the noise related to the variance of the phonon 
occupation. We provide an intuitive physical interpre- 
tation of this correction, which at large voltage agrees 
qualitatively with the asymptotic result of Ref. [l^ i.e., 
S ^ V*, and show that it generally leads to a significant 
increase of the noise above the phonon emission threshold 
in case of negligible phonon damping. Our predictions 
are relevant for currently ongoing cxperimentsj^ 

Model & methods. We consider the generic Anderson- 
Holstcin type of Hamiltonian for inelastic transport 
through a nanojunction H = He + u + Ht with 
the central part He = Ho + Hp^ + Heph, = 
J:^,Jh^J4dJ, Hpi, = nb^b,H,pi, = Y.^,,M^3dU3{b^ + 
b) and noninteracting electronic leads described by 
Hl,r = J2k,a=L,R ^o',kcl.k^a,k, tunncl-coupled by Ht = 

Ei,k,a=L,RiK,k^l,A + f^-C-) to the central part ("dot"). 
The states of the leads arc occupied according to the 
Fermi distribution /^(e) = (e/^^"^"^") + 1)"^ with /3 = 
l/ksT and /Iq, being the chemical potential of lead a. 
The applied bias voltage is — t^R = ■ We consider 
here only a single vibronic mode ( "phonon" ) with energy 
— other phonon modes simply contribute additively to 
the noise correction calculated below unless the phonon 



2 



resonances overlap, in which case a matrix generalization 
of our approach would be needed. 

To evaluate the noise, we employ the generalized 
Keldysh GF technique The key idea is to determine 
the GF of the dot G'^{t,t') = -ih~^{Tcdi{t)dj{t'y)x 
in the presence of a counting field A, introduced as a 
time-dependent fictitious parameter in one of the tun- 
neling matrix elements,— e.g., VI f. ^ VI j,e~^*^*^/^ and 
Vl*^ Vl*^e^^*^^^, with X{t) = ±A e R on the forward 
(labeled by ' — ')/backward ('-I-') branch of the Keldysh 
contour 1^ Knowing , the current and the noise (and 
in principle all the other current cumulants) can be di- 
rectly evaluated as described in detail in Refs. fill and [21I. 
In previous studies j^^—iii coupling to an external heat 
bath bringing phonons to equilibrium has been implic- 
itly assumed. Here, we focus instead on the limit of no 
external phonon damping, where all the thermalization 



comes exclusively from the electronic degrees of freedom 
and strong nonequilibrium phonon heating effects are ex- 
pected. 

Generalized kinetic approximation. To rigorously ad- 
dress the problem of nonequilibrium phonon heating, and 
its consequences for the transport properties of the junc- 
tion, one needs to include the influence of the electrons 
on the phonons by dressing the phononic GF with the 
polarization operator tlx. To the lowest order in e-ph 
coupling, this is given by the polarization bubbl o^^i^^'^^ 

nr(£) = -* / ^Tr[Mgr(£ + £')Mgr (£')], (1) 

with g^'' (cr, p = =f) being the Keldysh components of the 
dot's GF in the presence of the leads and of the counting 
field, but without e-ph coupling >iii2i 



Sxie) = 



el-h-^E„=L,Kr„[/a(e)-l/2] 
-iTLe-^[l - hie)] - iTr[1 - /^(e)] 



iTLe^fL{e)+irRfR{e) 
-£l + h-zEc.=L,flr„[/„(£)-l/2] 



(2) 



Here, the sign "stands for 2 x 2 matrices in Keldysh space and boldface letters indicate matrices in the dot electronic 
space indexed by the single-particle labels i and j, e.g., = {F^}, with F^(e) = 2tt J2k k^a k^i^ ~ £k,a) being 
the broadening due to coupling to lead a. Note that H^ is explicitly A-dependent and so is consequently the dressed 
phonon GF D\ given by the Dyson equation: 



Dxie) = 2n 
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with Vx{e) = (e^-n^ 
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In the limit of weak e-ph coupling l/Vx^s) shows 
two resonances around e ^ ±fl. Close to these res- 
onances, T^xie) is given by Vxie ^ ±fi) w 4f22 1^(£ ip 
O ± Aa(±0))2 + CA(±f^)/4], with Axis) = [H++(£) - 
n^-(£)]/2 and Ue) = 4H-+(e)n+- (e) - (H++(e) + 
H;^ (e))^- Neglecting the small real frequency shift 
A>(e) and approximating in the distributional sense 

s4t%4 ^ + 27r^(5(a;), where V stands for the 
principal value and S{x) is the Dirac (5-function, we arrive 
at the kinetic-limit expression - i.e., with zero phonon 
line-width - for the dressed phonon GF to the lowest 
(zeroth) order in e-ph coupling M: 
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where 



^2niY^ S{e + sn)J<fxie), 
(4) 
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nA+(e) 



= [Afxie)]'^ , is the A-dependcnt general- 



ized phonon 
to J^x=o{£) 
iU 



,+o(e)/2|ImH 



occupation. At A 

"W(e) + l/2 N{e) 

«(£)!, with H«(e) EE n^ro(£) -n^+o(£). 



0, it reduces 
where N{e) = 



is the nonequilibrium phonon occupation number, in 
agreement with Rcf. [l^- In the kinetic limit, N = N{n 
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coincides with the phenomenological results of Refs 
andlisl. It can be shown analytically that approximation 
(|4l) preserves charge conservation for all cumulants to the 
lowest (second) order in the e-ph coupling. 

Correction to noise due to nonequilibrium phonons. 
The effects of phonon heating on the current noise S can 
be evaluated within the formalism of Refjii, by replacing 
the free-phonon GF considered there with the dressed 
one, Eq. Doing so, we can split the correction to the 
noise at the lowest order in the e-ph coupling in two parts 
-I- 5ba- Here, S'av is the contribution due to 
coupling to a phonon with fixed average occupation, and 
it is simply given by the result of Ref . [ill with the thermal 
phonon occupation replaced by the nonequilibrium one, 
UBi^) N. In case of a single-level junction, it cor- 
responds to the phenomenological result of Ref. |^. The 
additional term S'ba represents the backaction of current- 
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driven nonequilibrium fluctuations of tfie pliononic occu- 
pation on tlie current noise itself and it reads 



5, 



ba 



with 



xi^)=^ I ^^^iP^s7ie-e')M. (6) 



After some algebra, S'ba can be rewritten as 



n' 



n' 



n' 



n' 



X {n{n + i)(n;+ + n'__ - 
+ (7v + i)(nV_-n'_+)-i(n; 



with u'^^ EE ^ul^in) 



. „. This expression is valid for 
an arbitrary junction with weak e-ph coupling and rep- 
resents one of the main results of this paper 

Further progress can be made within the extended 
wide-band approximation (eWBA)^ii^ where an ex- 
plicit expression for S'ba as a function of voltage and 
temperature can be derived. The full multilevel result 
is given in the supplemental material;^ while here we 
discuss for simplicity only the paradigmatic case of a 
single-level junction symmetrically coupled to leads at 
T = OfiiS In this case we obtain the following con- 
stituents of Eq. ^: N = e{\eV\ - Q){\eV\/n - l)/4, 

nv+ + n'__ - - n'_+ = i7eph(3 - 4r)ei//^, + 
n'_+ = -i7eph(i - 2r)eV/7T, n;_ - n'_+ = -i7eph(i - 

2nsign(l^)min(|ey|,r2)/7r, and |Imn-«(f7)| = -feph^/7r, 
with 7eph = M^7~^ /r^ being the dimcnsionlcss e-ph cou- 
pling constant and T € [0, 1] the clastic transmission 
coefficient. 

Analysis and discussion. We note first that S'ba is a 
strictly nonequilibrium correction, i.e., it is zero at = 0. 
This is consistent with the fluctuation-dissipation theo- 
rem, since there is no correction related to current-driven 
phonon fluctuations in the linear conductance. In the op- 
posite limit of large bias voltage, Sba is dominated by the 
term proportional to N{N + 1). Indeed this contribution 
grows like and for eV > fl, while the remaining 
terms of Eq. ([7]), as well as the "static" term 

Sav ^ in- 
crease at maximum as V^. The large- voltage behavior of 
Seph is then given by the leading terms of Sba 



Seph(ey » n) 



-{3~Ar)'^{eVfN{N + l). (8) 



This result can be fully understood in terms of a semi- 
classical mechanism related to slow fluctuations of the 
phonon occupation. The fluctuating occupation num- 
ber N{t) of a weakly driven oscillator is described by 
the master equation dPn{t)/dt = 7^[(n -|- l)P„+i(t) — 
nPn{t)] + 7^[nP„_i(t) - (n -t- 1)P„(0] for the probabili- 
ties Pn{t) that N{t) attains the value n at time The 



nonequilibrium rates can be microscopically evaluated 
as 7^ = in^^Q(r2)/fi,, 7^ = in^~p(r2)/fi, by comparing 
the rate equation d{N{t))/dt -(7^ - ^^){N{t)) + 7-^ 
for the mean occupation {N{t)) = T^^n nPn if) with 
the phonon-energy balance equation of Refs. |l2| and 
I15I . The stationary state has a geometric distribu- 
tion Pn{t 00) oc (7-t-/7|)" with the correct asymp- 
totic mean {N{t 00)) = (7;/7t - 1)"^ = N 
and the exponentially decaying connected correlation 
function {{N{t)N{Q))) ee {N{t)N{0)) - {N{t)){N{0)) = 



with the inverse relaxation time t, = 



7^ — 7^ = 2|Imn^(r2)|/S = 27ephri/fi.7r and variance 
{{N^)) = N{N + 1). At large voltage the correction 
to the mean current due to e-ph coupling is given by 
/eph(eF » n) « ^{N + i)(3 - 4T)eV ee i[N]2 
Since the relaxation time Tj-oi ^ ?i/(7ephf^) 3> h/D, ^ 
h/eV is the longest time scale in the problem, the noise 
at large voltage can be estimated assuming the cur- 
rent to follow adiabatically the phonon occupation N{t), 
Seph(eV^ » n) = dmNitMNiO)])) = {^f (3- 
ATf{eVf j^^dt{{N[t)N{Q))). Putting everything to- 
gether, we arrive at the microscopic result of Eq. 

The same line of reasoning can be repeated for a gen- 
eral multilevel junction without the need of eWBA. We 
can thus generally conclude that at large voltage (i) the 
correction to noise due to e-ph coupling is dominated 
by the backaction term Sba and (ii) this term is directly 
related to the diffusion of the energy stored in the os- 
cillator, which fluctuates slowly because of the random- 
ness of the driving tunneling events. The typical fre- 
quency of these fluctuation t"^^ is well above standard 
schemes for filtering out l//-componcnts (t~j^ ^ lOOGHz 
for h£l lOmeV, 7eph ~ 0.01), so that Sba will play a rel- 
evant role in noise measurements in systems with weak 
phonon damping. 

The backaction term Sba has also important conse- 
quences around the phonon emission threshold eV = fl. 
In particular, it affects the inelastic noise signal AS' = 
w\eV=Q+ ~ §v\eV=Q- ' '^hich is a quantity of prime ex- 
perimental interest defining the behavior of noise around 
the phonon emission threshold at low temperature. As an 
example, in the considered case of symmetrically coupled 
single-level junctions (in which AS' can be simply ex- 
pressed in terms of the transmission of the junction) the 
full result for noncquilibrated phonons Sav + Sba leads 
to AS' oc 247^ - 30r + 17/2, which changes sign at 
transmissions 7i,2 = 0.434 and 0.816. This should be 
contrasted with the values T1.2 = 0.146 and 0.854^ valid 
for completely thermalized phononi^ 

Finally, we have applied the eWBA multilevel results 
given in the supplemental material to the case of an 
atomic gold wire with unitary transmission and no ex- 
ternal phonon damping;^ as shown in Fig. [TJ Remark- 
ably, including Sba more than doubles the noise above the 
phonon emission threshold as compared to the "static" 
contribution Sav alone. For comparison we also plotted 
the noise in the case of fully thermalized phonons. An 
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FIG. 1. (Color online) Zero-temperature current noise in a 
perfectly transmitting atomic gold wire in the presence of an 
undamped alternating-bond-length phonon mode.— In this 
case Sciastic = 0. Solid red line, S = Sav + Shu', dashed blue 
line, 5av; dotted black line, case of thermally equilibrated 
phonons; inset, voltage derivative of the noise. 



intermediate level of equilibration simply interpolates be- 



tween the two limitsj^ This example clearly shows the 
need for taking into account the phonon backaction cor- 
rection in order to properly describe the noise above the 
phonon emission threshold. 

In conclusion, we have evaluated the full lowest-order 
expression for the current noise in the presence of an ex- 
ternally undamped phonon mode, identifying a purely 
nonequilibrium contribution due to the backaction of the 
heated phonon. At high voltage, this contribution dom- 
inates the noise and it can be fully understood in terms 
of slow fluctuations of the phonon occupation. Impor- 
tantly, this nonequilibrium correction also modifies the 
quantum behavior of the noise around the phonon emis- 
sion threshold and has, therefore, direct implications for 
ongoing experiments. 
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tion of the generalized kinetic approximation by avoiding 
unnecessary complex factors. 

The derivation leading to Eq. ([7| can be easily generalized 
to account for partial phonon equilibration due to a weakly 
coupled external heat bath, provided the appropriate A- 
independent phonon self-energy term is added to Eq. ([l]). 
For details, see the supplement.— 
See enclosed Supplemental Material. 
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Phonon heating also affects the crossover between enhance- 
ment and reduction of the conductance due to e-ph cou- 
pling. For the single level model this crossover shifts from 
T — 1/2 toT = 5/8 in the two opposite limits of strong 
and negligible external phonon damping. 



